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Abstract 

Limits at 90% c.l. on the square of the mixing strength |C/r4p between and 
a mostly isosinglet heavy neutrino with mass in the range 10-290 MeV/c^ are re- 
ported. The results were derived using the negative result of a search for neutral 
particles decaying into two electrons conducted by the CHARM collaboration in a 
neutrino beam dump experiment. Upper limits = 10""^ were obtained for neutrino 
masses larger than 160 MeV/c^. 
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Neutrinos may have Dirac or Majorana masses. In general the mass eigenstates 
(z^i, z/2, z^3, 1^4, ...) do not coincide with the weak (flavour) {ue^v^^VT^i-'s^ •••) eigen- 
states, but rather with a linear combination of them 

vi = ^UiiVi (/ = e, //, r, s, i = 1,2,3,4, ...) . (1) 

i 

Such a mixing could result in neutrino oscillations when the mass differences are 
small, and in neutrino decays when the mass differences are large. 

In this paper we report limits on the square of the mixing strength \ Uri\ between 
Vr and a heavy neutrino, 1/4, mostly isosinglet under the Standard SU{2)l gauge 
group, with mass in the range 10-290 MeV/c^. The limits were obtained using the 
negative result of a search for events produced by the decay of neutral particles into 
two electrons performed by the CHARM Collaboration in a neutrino beam dump 
experiment [1-4]. The decays of the neutral particles, produced in the dumping of 
400 GeV protons in a Cu target, were looked for in a volume located at a distance 
of -L = 480 m from the beam dump. 

The decay detector, shown in FigJT], has already been described elsewhere [5]. 
It had an empty decay region of D = 35 m length and 3 m x 3 m surface area 
defined by a veto scintillator plane (SCI) and a scintillator hodoscope (SC2). The 
volume was subdivided into three regions using two sets of four proportional tube 
planes (PI and P2) [6]. One module of the CHARM fine-grain calorimeter [6] was 
displaced to the end of the decay region. In order to improve the resolution of the 
shower angle measurement and to reconstruct better the decay point, three sets 
of four proportional tube planes (P3, P4 and P5) were installed in front of the 
module. Lead converters of 0.5Xo each were placed in front of PI, P2, P4 and P5. 
The detector was parallel to the neutrino beam line at a mean distance of 5 m, 
corresponding to an angle with respect to the incident proton beam of 10 mrad, and 
covered a solid angle of 3.9 x 10~^ sr. The signature of the neutral particles decaying 
into two electrons would be events originating in the decay region at a small angle 
with respect to the neutrino beam axis with one or two separate electromagnetic 
showers. 

The detector was exposed to a neutrino flux produced by 1.7 x 10^^ protons 
on a solid copper target [7] and 0.7 x 10^® protons on a copper target laminated 
with an effective density of one-third of that of solid copper [7]. In the combined 
exposures, 21 000 events were collected satisfying the trigger requirements of no hit 
in the scintillator planes SCI and a hit in at least four out of the six scintillator 
planes of the calorimeter module. The events were further selected requiring that the 
transverse co-ordinates of the shower vertex lie in a square of 2.5 m x 2.5 m centred 
on the detector axis and that the electron energy, E^\, measured in the calorime- 
ter module, be larger than 2 GeV. The events recognised as cosmic rays were also 
rejected. The remaining sample of 7185 events is dominated by inelastic scatter- 
ing of electron- and muon- neutrinos and antineutrinos producing hadron showers. 
Compared with the decay of neutral particles into two electrons, such events have 
a broader reconstructed angular distribution because of the intrinsic resolution and 
leakage effects. The regularity of the development of electromagnetic showers was 
used to distinguish further between the signal and the background events. In par- 
ticular, the distribution of the deviation of the reconstructed direction of the (two) 
shower (s) from that of the incoming beam and the fraction of the energy detected 
by the proportional drift tubes of the calorimeter module outside a narrow cone 
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Figure 1: Layout of the decay detector. 

around the shower were evaluated for the decay events by a Monte Carlo method 
[4]. No event compatible with the features of the decay of a neutral particle into 
two electrons was found. 





b) 



Figure 2: Feynman diagrams illustrating: (a) z/4 production from D^] (b) leptonic decay 
of T~ into a heavy neutrino; (c) decay of an isosinglet neutrino according to mode (|^). 



In this analysis we assumed that the 1/4 couples to matter essentially via mixing 
with Vt (|C/e4p , |C//^4p << |CA-4|^) and we studied the proton production of D^'s in 
the target and their decay into r. The heavy neutrino is produced in the charged 
current decays 

^ z^4 + r (2) 

and 

r ^ 1/4 + (3) 

The Feynman diagrams of the process (|2|) and of the decay t~ — > z^4 + /~ + vi 
{I = e, fi) are shown in Fig.s 2a and 2b respectively. On the basis of the assumptions 
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made, the isosinglet heavy neutrino decays only via neutral current interactions 
according to the modes 

1/4 ^ i/^ + e+ + e" (4) 

1/4 ^ Ut- + 1^1 + i>i (/ = e, /i and r) (5) 

^ + + fi~ (6) 

z/4 ^ z^T + 7r° . (7) 

The Feynman diagram illustrating the signal decay channel @) is shown in Fig. 2c. 
The channels contribute to the beam attenuation. The branching ratio of 

mode (^) is negligible and the decay (Q) opens for neutrino masses larger than the 
TT^ mass. The total decay width is then given by 



"tot = T {1^4 ^ I'r + 1^1 + I'l) + T {1^4 ^ Vr + + e ) + 



(8) 



For heavy neutrinos with mass larger than 290 MeV/c^ other decay modes open. 
The leptonic partial width is predicted to be [8]: 



F (1/4 Vr + vi + vi) +T {U4 Ur + + e~ 
{l + ~gl+gl) G2m5jC/,4|2 (1 - |C/,4|2 

1927r3 



(9) 



K 



where = 5l — 1 = —1/2 + sin^ and i^r, = sin^ ^'w is the weak angle. In 
this study the neutrinos were assumed to have Dirac masses and then K = 1. For 
Majorana neutrinos K is equal to 2. The leptonic partial width is dominated by the 
mode (^): 



F (1/4 ^ i/t- + e+ + e ) 



F (1/4 Ur + 1^1 + vi) + T {u4 Ur + e+ +e- 
The partial width for the decay (0) is predicted to be [9] 



0.14 



(10) 
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(11) 



For a given heavy neutrino mass rriy^ the number of the decay events (Q) expected 
in the detector is 



N 



where $ [Ey^) is the differential flux of heavy neutrinos, Pj 



(12) 

{Ei,^) is the 

probability for a heavy neutrino of energy Ey^ to decay in the decay fiducial volume 
according to reaction (^) , and e {rriy^ ) is the efficiency of the selection criteria of 
one or two electrons in the calorimeter based on the regularity of the development 
of electromagnetic showers and the collinearity between the (two) shower(s) and 
neutrino direction. The flux $ (E^^) is given by 



<^{E, 



U4J 



Nr 



BR (Ds ^U4 + t) A^l (Ey,) + 

+BR {D, -^Ur + T) BR {t^V4 + ...) Al^^l^ (E,,) 



(13) 
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The number of protons on the target corrected for the detector dead time, 13.6% 
for the sohd target and 21.7% for the laminated one, is A^p = 2.0x10^^. The fraction 
of proton inelastic interactions leading to a charged Ds is given by [10] 

^ _ AcA<r(D.)/.m.r-' ^ 2.98 . 10- . (14) 



Cinel CT; 



inel 



where the copper mass number is = 63.55. Linear A dependence is assumed for 
charm production. The used values of the inelastic proton cross-section, [11], 
of the ratio of the production cross-section for over the production cross-section 
for + D^, a {Ds)/(y [D] [12], and of the inclusive cross-sections for the produc- 
tion of D mesons, 0-2^=1*^°'^ [13], are reported in Table 1. The ratio a (Dg) /a {D) 
was obtained by the Beatrice experiment studying charmed particles produced by 
7r~'s of 350 GeV/c [12]. It is compatible with the results obtained by e'''e~ exper- 
iments at center of mass energies equal to 10 GeV and at mass [14]. The value 
of £7™!^!'=^°° was obtained by the NA27 Collaboration studying the production of D's 
in the interactions of 400 GeV protons in an H2 target [13]. 



Parameters 


Values 


Systematicerrors [%] 


^l^nl [mb] [11] 


769 ± 23 


3.0 


aiDs)/aiD) [12] 


0.12 ±0.03 


25.0 


^nuclcon ^] [^3] 


30.1 ±3.1 


10.3 


Br {D, ^Vr + r) [15] 


0.07 ±0.04 


57.1 


n [13] 


4.9 ±0.5 


4.0 


6 [13] 


1.0 ±0.1 


5.0 


Spectra of z/4 produced in r decay 


see text 


5.0 



Table 1: Values of the parameters used in the analysis and their contribution in percentage 
to the systematic error on the expected number of decay events (|^). The uncertainty due 
to the knowledge of the spectra of z/4 produced in r decay refers to heavy neutrinos from 
reaction (|^). 



In equation (|T3|) one has 

BR{Ds^VA + T)=BR{Ds^Vr + T)pDs\Uri? ■ (15) 

The value of the branching ratio of the Ds decay into a zero mass neutrino, 
BR (D^ Vt + t), is reported in Table 1 [15]. Its uncertainty dominates the sys- 
tematic error of this study. The factor pr,^ describes phase space and helicity effects 
[16]. In the case of r decay into 

BR(r^Z.4±...) = l^r4|'5]BR(T^Z/,+Xi)pt=Pr|t/r4|' , (16) 

i 

where BR [t ^ Ur + Xi) is the branching ratio of the r decay into a zero mass 
neutrino according to the considered mode i [15], see Table 2, and the p\ are factors 
depending on heavy neutrino mass. For i = 1-3 they take into account phase space 
and helicity effects [17]. For the modes 4-11 the p\.^s were computed using only 
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i 


Mode 


Branching ratio [%] 


1. 


T ^ ^1 + Uf, + Ut 


17.37 


2. 


T ^ e + Ue + Ur 


17.83 


3. 


T —>■ 7r~ + Ut 


11.09 


4. 




25.40 


5. 


r ^ 7r~ + X° + i/^ 


1.06 


6. 


r — >• TT + ztt" + 


9.13 


7. 


r ^ TT"*" + 27r~ + Ut- 


9.49 


8. 


r — >• 7r~ + 37r'' + Ur 


1.21 


9. 


T 27r- + 7r+ + 7r° + 


4.32 


10. 


r ^ 27r~ + 7r+ + X" + i/,- 


1.35 


11. 


+ 7 + + 


1.75 



Table 2: The r decay modes and corresponding branching ratio values used in the analy- 
sis 



phase space [18] . The numerical values of p'^ as a function of the neutrino mass are 
smaller than the ones of Ref. [19] by about 10%. The latter were obtained taking 
into account also helicity effects and the experimental width of the vector meson p. 

For a given neutrino mass, the acceptances of the heavy neutrino flux coming 
from the Dg [t], A^^ [^14]' ^"^^ corresponding energy spectrum normalized 
to one, (t)y^ (-£-1/4) [<Pl,^ (-E'1'4)] were obtained using a Monte Carlo simulation. The 
production of strange charm by protons was parametrized using the semi-empirical 
expression 

f{xF)^{l-\xF\re-^^^' , (17) 

where xp is the meson longitudinal momentum in the collision center of mass frame 
divided by its maximum value \/s/2, and pT is the meson transverse momentum. 
Since there are few experimental results available on the production of Df, the 
values of n and b were inferred from the measurements of D production. Assuming 
the hadronization process to be independent of the cc production mechanism, the 
parameters n and b are independent of the meson produced. Most measurements 
agree with a value of b equal to 1 (GeV/c)^^. The used values reported in Table 
1 were obtained by the NA27 Collaboration studying the production of D's in the 
interactions of 400 GeV protons in an H2 target [13]. Cascade production was 
neglected. 

The energy spectrum of heavy neutrinos from r decay is given by 

€AE..) = '-^, , (18) 

i=l 

where (f)^^ {Evi) is the normalized energy distribution of neutrinos produced in the 
decay mode i (see Table 2). In the case of leptonic channels, i = l and 2, the spectra 
were obtained using the matrix element 

|^|'~ (Pr-P.j (Prp.4) • (19) 



6 



The quantities p^, p^j^ , pi and p^_^ are the four-momenta of r, hght neutrino, electron 
or muon, and heavy neutrino, respectively. The spectrum of heavy neutrino pro- 
duced by channel 11 was obtained using phase space. The multi-pion decay modes 
were simulated using two models. In model (a) the spectra of channels 4-10 and 
their relative contributions as a function of the heavy neutrino mass were computed 
using phase space [18]. In model (b) channel 4 was assumed to be produced through 
the resonance p and channels 5-10 through the resonance ai. The resonances were 
assumed to have zero width. The p^'s in (|l^) are given by: 



Pr 



[l-yY + x{l + y-2x) 



1 -|- X — 



A 



2 + 2x-y 
{l-xf , 



(20) 



where x 



m„ / m,; 



for i = 4, {nip = 770 MeV/c^ 



m 



^^/m^ for i = 5-10 (m^j 



= 1260 MeV/c^) and y = Jm^. The values of Table 3 were computed using the 
average of the spectra obtained in the two models. The systematic error in Table 1 
reflects the differences of the spectra. 

The decay (H) was simulated using the matrix element [8] 



- gl {PvA • Pe- ) {Pu^ ■Pe+)+ 9r {PvA ■ Pe+ ) {Pi^r ' Pe- ) 



(21) 



which neglects the electron mass. The quantities p^- , p^+ and p^^ are the four- 
momenta of electron, positron and tau neutrino, respectively. In the center of mass 
of the decaying heavy neutrino the four- vector p^^ is given by [20] 



Pui = {m. 



va 1 



-m 



U4, 



\h\fi) 



(22) 



where r/ is a unit vector parallel to the direction of the heavy neutrino in the rest 
frame of the particle decaying into neutrino, Dg, or r, and \h\ is the absolute value 
of the neutrino (antineutrino) helicity. In the case of heavy neutrinos coming from 
Dg decay, the values of \h\ obtained in Ref. [15] were used. As the polarization of 
r produced in Dg decay is negligible, = was assumed for the heavy neutrinos 
produced in r decay. The acceptances and the mean momenta of decaying heavy 
neutrinos expected to be detected in the detector are reported in Table 3 for different 
values of neutrino mass. The efficiency of the cut > 2 GeV is about 85% for 
heavy neutrinos coming from Dg and larger than 95% for the ones coming from r. 
In Eq. (12) the probability for a heavy neutrino of energy Eu^ to decay in the 



decay fiducial volume is given by 



r(i^4 



+ e'^ + e 



r 



tot 



where A = (7/5c)/rtot is the heavy neutrino mean decay path (7 
Py^l Ey^). According to Eqs. (^ and ([Tl|) A depends on |C/t-4|^ [l - 



E,, 



'-'vaI ^i^a 



(23) 



The quantity e {m^^) in Eq. (12) is the efficiency of the selection criteria based 
on the regularity of the development of electromagnetic showers and the collinearity 
between the (two) shower (s) and the neutrino direction. The values, obtained using 
a detailed Monte Carlo simulation of the detector response, decrease with increasing 
heavy neutrino mass. It ranges from 91% for rrip^ = WMeV /c^ to 65% for rrip^ = 
290MeV/c2 [4]. 

Since no decay event was detected, upper limits at 90% confidence level on 
|f7-r4| were obtained in the neutrino mass range 10-290 MeV/c^. The hmit value 
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r T\ /T \T 1 21 

rriv^ [MeV/c J 


A^;^" X 10 


{Pu:} [GeV] 




{PQ [GeV] 


10 


3.39 


14.39 


3.49 


47.14 


50 


3.50 


14.14 


3.54 


47.13 


100 


4.11 


12.90 


3.61 


46.86 


150 


6.00 


10.39 


3.73 


46.48 


190 


8.00 


8.38 


3.87 


46.44 


250 






4.04 


46.03 


290 






4.24 


45.62 



Table 3: Acceptances and mean momenta of decaying heavy neutrinos with Eei > 2 GeV 
for different values of the neutrino mass. 



of N = Ni = 6.42 events was used in Eq. (12). Since the contribution to the 



systematic error of the uncertainty on the spectra of heavy neutrinos coming from 
r decay is neghgible, Ni does not depend on the neutrino mass. It corresponds to 
an average probability of observing no events, {Pq (Ni)), equal to 10% [21]: 

+00 

{Po m) = J fi^-^ ^0 ^ (^/; ^) = 0.1 , (24) 

—00 

where / (0; A'^;') = e~^i is the Poisson probability of obtaining zero events. In 
the case of negative A'^;', / (0; N[) = 1 was used. The probability density function 
W {N'l^; Ni,a) takes into accounts the systematic errors summarized in Table 1 and 
was assumed to be a Gaussian distribution. Combining in quadrature the uncer- 
tainties reported in the table, one gets a/Ni = 0.64. In the case of no uncertainty, 
W {N[;Ni,a) = 6 {N^ - Ni) and the integral gives Ni = 2.30: the upper limit at 
90% confidence level for the mean value of a Poisson distribution in the case of zero 
observations. The chosen value A'^; = 6.42 is a safely conservative number: taking 
for instance W {Nl; Ni, a) equal to a log- normal distribution to avoid negative val- 
ues of Nl [21], one obtains A'^; = 3.21, a factor 2 smaller than the value used in the 
analysis. 

The upper limits obtained at 90% confidence level on |C/t-4|^ values, as a func- 
tion of shown in FigJ^, together with previous results [9, 22]. Limits 
on |?7t-4|^ were also obtained for neutrino masses larger than 140MeV/c^ assum- 

1 2 1 2 1 2 

ing \Ue4\ = \U^4\ = \Uri\ from the upper bounds on the rates of the decays 
T~ (a*^) vr^vr" [23]. Limits on |C/e4|^ and |?7^4|^ are reported in Ref. [15]. 

In conclusion, the negative results of a search of decays of neutral particles into 
two electrons performed by the CHARM Collaboration in a neutrino beam dump 
experiment, allowed limits to be set at 90% c.l. on the square of the mixing strength, 
\Uri\ , between Vj- and a mostly isosinglet fourth neutrino, 1^4, having a mass in 
the range 10-290 MeV/c^. Values of ^10"'' were obtained for masses larger than 
160MeV/c2. 
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Dolgov, S.N. Gninenko, P. Loverre, M. Mangano, J. Swain and K. Winter for many 
interesting discussions and comments. We would also like to thank S.H. Hansen for 
advice on the astrophysical and cosmological limits. 
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Figure 3: Limits at 90% confidence level on the square of the mixing strength \Ur4\ of 
the r-neutrino with a fourth neutrino mass eigenstate mostly isosinglet : (a) upper limits 
from this study; (b) the NOMAD upper limits [22]. The lower limits (c) from SN1987a 
and (d) from the Big Bang NucleosjTithesis constraint AN^^^ < 0.2 are reproduced from 
Ref. [9]. Upper limits of {U^if < 10"^ from SN1987a and of < 10~^° - 10'^^ from 

Big Bang Nucleosynthesis have also been derived for the corresponding mass range shown, 
respectively [9]. All the limits were obtained assuming that neutrinos have Dirac masses 
{K=l). 
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